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I. INTRODUCTION 


This study was part of an ongoing project to produce an 
empirical heat budget parameterization for utilization in the 
Fleet Numerical Weather Central (FNWC) primitive equation 
prediction model. The project will eventually examine the 
heat balance of the earth-ocean system using representative 
mid-seasonal data from each of the four seasons of the year 
1973-74. This particular study covers 15-16 October 1973, 
using gridded FNWC temperature and humidity fields, at con- 
stant pressure levels, for this period. 

The parameters which influence the radiative calculations 
most strongly are the two cloud layers of fractional amounts 
CL(1) and CL(2), at 600 and 900 mb, respectively, at each 
grid point. The fractional amounts are specified using the 
large-scale formulations developed by Smagorinsky (1960) as 
also was done by Jenks (1974). The possibility of two cloud 
layers introduces substantial complexity to the radiative 
physics compared to the one-layer cloud model previously de- 
Seribed by Martin (1972). The adaptation of the radiational 
model to FNWC gridpoint soundings is based on an unpublished 
fmanuwseripty (Martin, 1974). The Martin manuscript in turn 
has similarities to the multi-cloud heating packages of 
Arakawa (1972) and Rodgers (1967). 

Heat-balance computations were made at the ocean surface, 


at the tropopause, and for the atmospheric layer surface to 


16 





200 mb. A major purpose was to determine whether the cloud- 
cover parameterization previously adopted by Jenks (1974) 
lends itself to realistic budgetary computations (as compared 
to the results of others) for Oct. 15-16 over the Northern 
Hemispheric oceans. In studying the net heating effects at 
the ocean surface, the evaporative and sensible-heat trans- 
fers were modeled after Kesel and Winninghoff (1972), and 
Kaitala (1974), but with some modifications to be described 
aber (Section V). 

To assess the effects of the cloud cover parameterization 
on the heat transfers of this study, gridpoint data along 
four meridians were selected from the FNWC analyses of 15-16 
October 1973. Three of the meridians are from the Pacific 
Ocean area and the fourth is from the Atlantic Ocean (see 
Fig. 2). Most of the gridpoint soundings sampled were in 
the Northern Hemisphere. 

The five sigma (o = p/1) analysis levels of the FNWC 
primitive equation model were adapted for data input in the 
vertical scale of this study. The k-levels (k=1,...,10), 
where k=100, are used for the purpose of simplified refer- 
encing of the data levels in the equations derived and are 
SHowiermm tie. ft. The cloud layers, when they exist, are con- 
Sidered to be in two layers, a high-level cloud layer from 
k=4 to k=6 and a low-level cloud deck from k=8 to k=9, with 
amounts CL(1) and CL(2), respectively. In pressure-scaling 
the water vapor and CO, absorber masses a realistic method, 


after. Smith (1966), was used to extrapolate water vapor 
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(8 10) om = I enaaieaadtestataeas ae 09 
| tt 
Figure 1. Five-layer radiative sounding used in this study. 
Levels are identified by their values on the k-Scale, while 
layers are identified by their level boundary indices in 
parentheses, c.g. (8,10). Pressure-scaled water vapor and 


CO. mass increments U and C, respectively, are integrated 
with respect to the surface and introduced at even levels 
Whtbew tie Mixing ratio, q, is formulated at odd bevels. The 
femperarture, J, 1S expressed for all levels. 
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mixing ratio values above the level of k=4, where radiosonde 
moisture data is normally not available. The net flux of 
terrestrial radiation computations make use of empirical ab- 
sorptivities of water vapor and CO, (Sasamori, 1968), which 
are incorporated in the Martin (1974) manuscript. The solar 
radiation package in Martin (1974) has been modeled after 
Arakawa (1972), except that only two cloud layers are con- 
sidered possible, as contrasted with three in the Arakawa 


version. 
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If. DATA PREPARATION 


The data used in this study was in the form of gridpoint 
soundings chosen along four meridians of FNWC 63-by-63 grid- 
mesh temperature and humidity data fields to insure minimum 
over land data. The three meridians (and their respective 
number of soundings) selected over the Pacific Ocean were 
tgeateaq at 125 W - (25 soundings), 170°W —- (25 soundings), 
and 145°E - (17 soundings). The Atlantic Ocean meridian was 
35°W - (26 soundings). Fig. 2 shows how the meridians were 
located on the FNWC polar stereographic grid. This procedure 
of choosing data lines along diagonals or along J=32 (through 
the pole-point) of the FNWC data-fields obviates the neces- 
sity of interpolating spatially to gridpoints. Data along 
line 3 in the Pacific were not extended southward of (9,55) 
because they fell over land masses (New etenes and Northern 
Australia) where the surface temperatures were too high to 
be representative of the oceanic values. 

The vertical soundings on these data lines were taken 
from the original FNWC 63-by-63 gridpoint analysis of tem- 
perature at the surface and at nine other standard levels up 
PomlOoommulltpars £Lor QOOOGMT, 16 October 1973 for the Pacific 
Ocean data cases and 12 hours earlier for the Atlantic Ocean 
data cases. In addition to these standard level values, a 
ie value was listed which is the FNWC calculated temperature 
(see Kesel and Winninghoff, 1972) for the top of the turbu- 


Jent layer of constant sensible heat flux. This level, 
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(1,63) 


t A~145°E 


(1,32) 


A= 170° W 





(1,1) _ bs 
i=25 i=38 ge (63,1) 


Figure 2. FNWC polar stereographic grid and meridians (lines 
2, 3, and 4) selected for study. The longitudes A are 


shown for each meridian as well as the extent considered of 
each merialan. 


au 











which iS approximately 20 meters above the surface, was des- 
ignated by a code label of 9999 under the column labeled 
pressure (Table I(a)). The six moisture values of this table 
were given by the vapor pressure at sea level, and by dew- 
point depression for levels from 925 mb to 400 mb. In addi- 
tion, the parameter e, was provided by FNWC as the vapor 
pressure at the top of the constant flux layer at the level 
indicated by the code 9999 (Table I(a)). Table I(a) gives 
an example of the original sounding for gridpoint (1,1) of 
the FNWC map for the given period. 

iMmondereGo periorm radiative calculations, it is neces- 
sary to have water vapor and CO, absorber masses and cloud 
amounts within the required k-level boundaries (Fig. 1). To 
do this, a systematic method was applied to arrive at the 
temperature and mixing ratio distributions (Table I(b)) as 
detailed in the following subsections. It should be noted 
that all soundings start at sea level, but as a reasonable 
approximation surface pressure was considered to be 1=1000 mb. 
Therefore the k-levels of Fig. 1 become 1000., 900., 8O00., 


., 200., 100., 0.0 mb, respectively. 
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ie sta. Bxample of an original FNWC gridpoint sounding 
fomepoilnt (1,1). The humidity parameters for the surface 
and the tope of the constant flux layer (9999) are vapor pres- 





sures and those between 925,...,400 mb are dew-point depres- 
sion. 
Pressure (mb) TEC) Humidity Parameters 
surface 20.9 25.6 mb 
925 19.3 Silo C 
850 Leo Sae € 
700 6.5 See 
200 -10.2 baie 
400 -21.9 Pi MAG: 
300 —-36.8 
Za0 -45.8 
200 -56.0 
150 -67.3 
100 -80.6 
9999%* T= 24.7 e. = 26.3 mb 
Table I(b). Example of the corresponding radiative sounding 


with temperature and mixing ratio listed at k-levels. 


* Values at top of constant flux layer as used 
ing package (Kesel and Winninghoff, 
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OG). 


Pressure (mb) T(°K) Mixinee Rat To. 2) ke) 

1000 298.7 15.92 
900 290.9 11.45 

800 Zoom 
700 20 aes, Ono 

600 272.4 
500 263.0 D3 
400 Zoe 1.42 
300 236.4 0.56 

200 ZA 2 2 
100 192.6 0.03 
0 LOZ. 6 0.00 


in FNWC heat- 











A. TEMPERATURE TERMS 

The temperatures were listed at each mandatory level be- 
tween the surface and 100 mb. The temperature T,9 was set 
equal to the listed surface air temperature. At the top of 
the atmosphere, the temperature was assumed to be isothermal 
from 100 mb to P=0 mb, with the mixing ratio approaching 
zero at P=0O mb. The radiative sounding temperatures for the 
remaining k-levels were obtained from either their corre- 
sponding listed temperature or by a three-point Lagrangian 


interpolation scheme of the form: 


Ti. = To (CPL Pi )(P,-P2 ))/CCPo-Pi ) (Po-Pi )) 
. Ti CCP) Po CP, -P2))/CCP1-Po )(Pi-P2 )) (Zari 


+ Te (CP, Po CP, -P1))/C(CP2-Po )(P2-Pi )) 


Here k is the pressure level desired for T, interpolation, 


k 
P, lies above Pio and P,, P, lie below level k. The param- 
eters T,, T,, and T, are temperatures at the same levels as 
P,, P; and P,, respectively. As an example, if k=9, P, and 
Po were taken as 925 mb and 1000 mb, reSpectively. 

The ue temperature Was MOL Ubalized in the interpolation 


scheme described above, but was carried along as a twelfth 


temperature for later use in the heat budget computations. 


B. MOISTURE TREATMENT IN THE SOUNDING 

The moisture parameters from the original FNWC soundings 
(Table I(a)) were converted into mixing ratios at each of 
the original sounding levels. At k=10 the surface vapor 


pressure, ore was used to calculate the mixing ratio as 
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follows: 
Gdiooo = wOZLOT e,yyp/ 1000 (2-2) 


where e,,,. is in millibars. To obtain a mixing ratio tor 
the top of the constant flux layer, Equation (2-3) was used 


to compute q, from e, 


q. = (e,/1000) .62197 (2-3) 


where cy 1s the vapor pressure at the top of the constant 
mie layer in mb. 

The mixing ratio at level P is calculated at the remain- 
ing original data levels by the use of Equation (2-4) which 


is taken from Fleagle and BuSinger (p. 63, 1963); 


A-B/T) -B(T-T_ ) 
621.97 e' 
q(P) = 6+ Eee (Zea) 
Pp hell 
D 

A = 21.656 

B = 5418.0°K 

T = temperature °K at level P 

Ty = dew point (°K) at level P 


The calculated q-values were then interpolated to the k-lev- 
els using an equation similar to Equation (2-1). An example 
of the interpolation results for gridpoint (1,1) is shown in 
habe. 1(b). 

Since accurate humidity information is not generally 
available for the upper troposphere, a logarithmic extrapola- 


Elon formula was used to find q-values at k=1, 2, and 3. 


a(P) . (Py) (265) 
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This formula as described by Jenks (1974) is similar to that 
developed by Smith (1966). It is solved, using a "best-fit" 
between variables y = log ae and x = log ant for the param- 
eter X. The six q-values of the original soundings corre- 
sponding to P = 1000, 925, 850, 700, 500, and 400 mb are 
used to determine i, which is then used in Bacon (2-5) to 
extrapolate the q values to P = 300, 200, and 100 mb. 

A correlation coefficient was statistically computed in 


connection with each A-profile of the form: 





The correlation coefficient enabled the monitoring of the 
usefulness of the A-profile technique and was generally 


found to be in the range of 0.95 to 0.99. 


C. PRESSURE-SCALED ABSORBER MASSES 

Having computed the mixing ratios at each odd k-level, 
it is now possible to compute the pressure-scaled water va- 
por absorber masses in each layer. In Figure 1, these water 
vapor layer-masses were denoted by U, however here we will 
use notation which is more representative of the particular 
layer under consideration. For example, the pressure-scaled 
water vapor mass in the lowest layer (8,10) is given by 


qgAP Pg aie 


0 
Ee) = see ‘j013. 25! 


(2-6) 
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where the Moller and Raschke (1964) pressure-scaling exponent 
0.72 has been employed. For the layer (6,8), an analog to 
Equation (2-6) was used by replacing the parameters with the 
subscript 9 with their counterpart parameters with the sub- 
script 7. A Similar procedure was followed for the remain- 
ing layers (4,6), (2,4), (0,2). In each of the five layers 
AP = 200 mb. 

The algorithm for computing the integrated water vapor 


mass above the earth's surface is as follows: 


WS) =O. 

Wes = UClO), + UCs, 10) 

U(6) = U(8) + U(6,8) 

U(4) = U(6) + U(4,6) (27) 
Weep = UK4)) + U(2,4) 

UC} = U(2) + 0.5U(0,2) 


WMGODee = UC2Z) + UlO,2) 


A similar computation scheme is used for the carbon di- 
oxide scaled mass. Here, however it is customary to repre- 
sent the pressure-scaled mass in terms of the normal 
temperature and pressure (N.T.P.) volume over a cm* of the 
earth's surface. Thus for the lower layer: 


AP Py eee 


0 
= se 7) tees ee ee 


(2513) 
The fractional amount 3.14 x 10-* is the constant mixing ra- 


tio of COs by volume, and = 





is the N.T.P. thickness of the 
S 
layer AP = P,o9-Ps. The quantity P. is constant density ina 


Standard homogeneous atmosphere of depth H where T=273.16°K 
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everywhere and P. = 1013.25 mb. The parameter H, also gen- 
erally known as the scale height of the homogeneous atmos- 


phere, may be further defined as: 


H oe = 7,995 x 105 cm 


Therefore Equation (2-8) may be rewritten as 


AE rE 


0. 
BEER Pat = 
H(7013.95) {013.95 (2-9) 


CeSe1o) = 3.14 x 107° 


All other layer thicknesses of CO) are obtained in a similar 
manner to Equation (2-9) by substituting the appropriate 
subscripted indices 7,5,...,1 on the pressure P. 

The algorithm for computing the integrated carbon diox- 


ide mass above the earth's surface is as follows: 


C(10) 0.0 

C(s) = C(10) + C(8,10) 

C(6) = C(8) + C(6,8) 

C(4) = C(6) + C(4,6) C20) 
C(2) = C(4) + C(2,4) 

Cees 7C(Z) + 0.5C(0,2) 


C(O) = C(2) + C(0,2) 


D. CLOUD PARAMETERIZATION 

In order to determine the fractional cloud cover amounts 
in the layers (4,6) and (8,9), the saturated vapor pressure 
must be calculated at levels k=5 and k=9, respectively by 


the following equations: 


(A-B/Ts5} 
O29 1 ee 
Ig (500) 500 ( a ) 
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Gi=B/ Ts) 
.62197 
q,(900) = ett (2-11b) 


where A = 21.656 and B = 5418.0°K as before in Eq. (2-4). 
Knowing the saturated vapor pressure at levels k=5 and 
k=9 enables the determination of the relative humidity at 


those levels: 


_. g(9500) 
Sagoo) _ 


The relative humidities were then substituted into 


omagorinsky's (1960) formulation for fractional cloud amounts: 


CL(1) Zor Che) ) =O. 7 (2-13a) 


iD 


Soom Gino) — 2.0 (2-13b) 


These formulations were employed by Jenks (1974) and are 
also to be tested in this study in which a different season- 
al sample iS considered. 

These fractional cloud amounts are a function of large 
scale influences only and do not consider small scale con- 
vective activity, seasonal conditions or latitude. The use 


of such clouds is for radiational study only. 
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III. TERRESTRIAL RADIATION 


A. THEORETICAL AND EMPIRICAL BASIS 

Empirical formulas were developed by Sasamori (1968) for 
computing emissivities for long-wave flux calculations asso- 
ciated with the NCAR General Circulation Model. These em- 
pirical emissivity formulas can be related to the values 
built into Yamamoto's (1952) well known Radiation Chart, 
which was developed from radiative transfer theory. This 
chart has proved to be very accurate for computational pur- 
poses, and was used by Jenks (1974) as a schematic guide for 
integration of the radiative transfers through the various 
layers of interest in the soundings discussed in Section II. 

The necessary emissivity formulas of Sasamori for utili- 
zation in the FNWC heating package are those which will yield 
expressions for: 


F,o* = Net IR flux at earth surface 


Fe* = Net IR flux at level k=6 
Fo* = Net IR flux at level k=2 
F610 = Net IR flux divergence in the layer (6,10) 
ieowe-elet IR flux divergence in the layer (2,6) 


A workable scheme for making theSe computations is presented 
Polewetorevarious combinations of cloud cover [CL(1), CL(2)]}. 
For example, in order to compute 


Bio 


Piot = Bio - a €,(U,C,7T)dB (3-1) 
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where B, is Che stetan-Bortzmamn blacksbody flwx at T one 


k? 
must have a good representation of the emissivity Cae as a 
function of both water vapor and carbon dioxide absorber 
masses (U and C, respectively) along the sounding which is 
miempath OL integration. For the integral evaluation of 
Eq. (3-1), Sasamori (1968) proposed Ea eraidt sandapenaene 


emissivity formulations ae as exemplified for the absorber 


masses in layer (8,10) as follows: 


E,- (8,10) = {.240 log, ,)[U(8,10)]+.622}+.07262{ (1. 
=, S25 56 PUCSmnO + O2e6)° *° )[ logy peCS510) 
+1.064) } (3-2) 


hornulas Of the type given by Eq. (3-2) apply for T > 210°K. 
The quantity in the first brace of Eq. (3-2) indicates the 
water vapor emisSSivity alone and that in the second brace is 
the additional contribution by that part of the CO, emissiv- 
ity not already included in the water vapor emissivity ex- 
pression. 

Eq. (3-2) was used in the calculation of the atmospheric 
Samcoivables Over all layers (6,10), (4,10),....,(1,10), for 
whach Cumulative mass values U(6,10), C(6,10), etc., were 
available from the radiative soundings. Thus Eq. (3-2) 
proved very versatile. For completion of the integration 
from level k=1 to the origin (B=0, ce=1.0) of Yamamoto's Ra- 
eeeation Chart, the UCT, ), c(T,) soundings were taken to be 
iaGtiermalwat T=T, to the top of atmosphere (to level "0"). 
This leaves only the temperature-dependent emissivity € 


WC 
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etgee tO), 1, )] which applied for T < 210°K, also cOmputed after 


Sasamori (1968) by means of 


= U(0,10),7T; g. 347, £353 log: oU(0,10)-. 44] 


wc 


emo ee log, ,»U(0,10)-. 705] 


x 


x 


[8.0/(.353 log, ,U(0,10)+3.56)] 


.07262{ (1.-.62556[U(0,10)+.0286}° 2°) 


om 


Ketel@Oeante (0.10 )+1, 064 |} (3-3) 


The first term is temperature-dependent, as well as dependent 
upon the summed absorber masses relative to the reference 
level (in this example, level k=10). The term in the brace 
involves essentially the same form of the CO2 correction to 
the total emissivity as seen in Eq. (3-2) except that it re- 
fers to the full layer (0,10). Eq. (3-3) is applicable with 
appropriate final parameters U(0O,k), C(0,k), and T, to the 
determination of the downward flux contributions through 
level k for all of the reference level cases to be consid- 


ered below, i.e. k=2,6,10. 


B. NET FLUX Fj, 0* 

The radiative soundings were evaluated in the form of 
the parameters U(k,10), C(k,10), and T,, level by level. 
Also listed for each radiative sounding are the cloud cover 
Popamewers CL(1) and CL(Z), which in general are non-zero. 


The grid area may be visualized as composed of areal frac- 


tions or weights defined as follows: 
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I 


W(0,0) Ci-CL(1))C1-Cil(2)) : totally clear fraction 


W(1,0) = (CL(1))(1-CL(2)) : overcast in upper layer only 

(3=4") 
Mneee ey) =) (1-CL(1))(CL(2)) : overcast in lower layer ondy 
WC1,1) = (CL(1))(CL(2)) : overcast in both layers 


The composite radiative effect of the true cloud cover 
over the grid area is then reconstructed by multiplying the 
results of the respective cases defined above by the appro- 
priate weight factors. The four cases will now be defined 
in order. 


1. Net Flux F,,* with Clear Skies 


Defining the Stefan-Boltzmann blackbody flux at Ty, 


as B,, where B, = 1.170403 x lO? a ly(day)~', the™net fix 
at the surface under the clear sky case is defined as fol- 


lows: 
Beno co, 0) = Byo-. Ste, (8,10) (Bi o-Be )+(Bs-Be le, (6,10) 
te (8,10) ) tle, (4,10) te (6,10) ] (Be-B, ) 
Coa) 
+(B,-B2 )[e, .(2,10)te, (4,10) ] tLe, (1,10) 


te, ,(2,10)](B2-Bi )+@,.[ (0,10), T1]Bi} 


2. Net Flux F,,* with Overcast Clouds in Upper Layer 
Only (k=4 to k=6) 


rio (tO) = (Bi o-Bs ){1.-.5fe, (8,10) (Bi o-Bs ) 


(3-6) 
+(e, (8,10)+e, (6,10) )(Be-Bs )1} 
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2.) Nelms Mig t with Overcastsin, Lowen, bay.cm Only 
(Kaeo k=9)) 


Brot (On k= (Bi o-Bs )[1.-.5€,.(9,10)] (3-7) 


4, Neto Plux F,,* with Overcast Clouds in Both bayers 
CIRO manG as. © ) 


Gs) = (By, 0-Bs )[1.-.5€,,(9,10)] (3-6) 


It should be noted that Eq. (3-7) and Eq. (3-8) are identi- 


cal. 


a Composite F,,* Calculations 


The composite results are obtained by multiplying 
the above defined four cases by their respective weights as 


defined by Eq. (3-4) to obtain 
Fio* = [1-CL(2)]{(Bio-Be )-.5€, (8,10) (Bi o-Be) 
Le, (8,10)+e, (6,10)] (Be-B)} 
+(1-CL(2))(1-CL(1)){Be-. 5[(e, (6,10) 
+e, ,(4,10))(Be-By +(e, (4,10) te, (2,10) ) (By Be ) 
+(e, (2,10)te, (1,0))(B2-Bi )+€,((0,10) ,T: Bi }} 
+CL(2){(Bio-Bs)[1-.5e,.(9,10)]} (3-9) 


It is clear that Eq. (3-9) reduces to F,,9*(0,0) when CL(1) = 
Ci(2) = 0, and to Fio*(1,1) when: GlG@ly = Clh(2) = 1.0. Ina 
similar manner the other two cloud cases may be recovered by 

setting the weights to CL(1) = 1.0, CL(2) = O and CL(2) = 1.0, 


CL(1) = 0, respectively. 
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The above described procedure for obtaining the net 
flux can Similarly be extended to levels k=6 and k=2 as 
listed in the following subsections. 

The composite F,)* values will be compared in Section 
IITI.F with climatology (Malkus, 1962) to ensure the reason- 


ableness of this calculation. 


Ce NED FLUX I, * 

The individual cases for F,* can be obtained from the 
composite F,* as described above for F,,)*. The numerical 
from using the Sasamori type emissivities [Eqs. (3-2) and 


(3-3)] for the composite F,* is as follows: 
Fe* = {1-CL(1)1{B,-.5le, (6,8)(Bs-Be te, (4,6) (Be-B, ) 
+(e, (4,6) te, (2,6) )(By-Bz +(e, (2,6)te, (1,6))(B2-Bi ) 
a eme@aOre) ty) Bu) + (1-CL(1) )(1-CLC2)) 1 (By 0-Be )[4. 
-.5(e, (6,8)te, (6,10) )]}+CL(1){ (BeBe )f1-.5(e, (6,8) )]} 


+CL(1)(1-CL(2)){(Bio-Bs )[1-. 5(e, (6,8) +te,,,(6,10))}} 
(3-110) 


D. NET FLUX F2* 

In a2 similar fashion to F,,* and F,* the net flux at 
Mevel Kk=2 can be calculated from the radiation chart using 
the trapezoidal rule as described by Jenks (1974). Again 
the numerical formulation of the individual cloud cases can 
be obtained from the composite net flux formula as described 


in Subsection III.B.5. The composite F2* formula is 





Fo* = [1-CL(1)]{Be-.5{c, (2,4) (By-B2 +(e, (2,4) 


Fe (2,6) )(Be-B, )+(e, (2, 6)te (2,8)) (BeBe) (3-11) 


+e ,(1,2) (Ba -Bi +8, ( (0,2), 7; )By]}+(1-CL(1))(1-CL(2)) 
x{ (Br 9-Bg )[1-.5(e, (2,8) te, (2,10)] }+CL(1){B, 
-. Sle, .(2,4) (By-By te, (1,2) (B2-B, )+€, ((0,2),T: Bi J] 


If the terms Eo (00,2),T1 By and Eo (1,2) (Bz -B, ) are set 
to zero in Eq. (3-11) a new parameter can be defined as the 
net outgoing long-wave radiation and labeled F,. This pa- 
rameter was computed to compare with satellite climatolog- 


ical observations over the oceans (Raschke, et al., 1973). 


E. APPLICATIONS TO HEAT BALANCE COMPUTATIONS 

For computation of the radiative "balance" at the ton of 
the earth-tropospheric system only F,* of Eq. (3-11) and the 
total insolation absorbed below k=2 are required. At the 
surface, a radiative balance is computed at each gridpoint 
using the absorbed insolation and the net flux F,)* of Eq. 
(3-9). Also at the surface a "convective" balance is com- 
pea by ancludinge the contributions of evaporation and seéen- 
Sible heat transport (see Sec. V). Net flux F,,9*, evaporation 
and sensible heat transport are all considered negative con- 
tributions at the surface. The atmospheric radiative bal- 
ance requires knowledge of the long-wave cooling effects 


caused by the flux divergences F26 and F610 as defined by: 
F26 = Fo*¥ - FF. * (3-12a) 


F610 = Fe* - Fi ot (3-12b) 
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As will be seen in a later section, these parameters act as 
destabilizing influences in their respective layers when as- 
sessed in relationship to other processes in cross sectional 


display form in Section VI. 


F. COMPARISONS WITH PUBLISHED RESULTS 

Table II gives a visual comparison of the outgoing long- 
wave radiation, F, (as contrasted with F,*), found by this 
study as compared to the mean values published by Raschke 
et al. (1973). Also displayed are this study's global net 
flux, F,.*, computed using Eq. (3-9) as compared to the an- 
nual mean deduced from Malkus (1962, Figs. 9,10). The 16 
October 1973 results have been interpolated, gridpoint by 
gridpoint, to 5° latitude increments for each of the four 
lines (see Fig. 2) and then latitudinally averaged to arrive 
at the "global" value for each latitude. The units of the 
values from Malkus @ctually due to Budyko (1956)) have been 


-! to those of cal cm-*min-! used 


converted from kcal cm-?yr 
in this study. The global mean values listed at the bottom 
of each column were derived by cosine weighting the latitu- 
dinal values listed in Table II (see also Sec. VII). 

It should be noted that the values of F,. and F,,.* found 
by this study for latitudes 20°S, 10°S, and 60°N are averages 
of only two "sounding" lines while the other latitude values 
are averaged over all four meridional lines, the difference 


being due to the number of soundings available near the ex- 


tremities of the FNWC cross-sections. 
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The consistently smaller values of net flux, Fi9*, found 
by this study might be explained by the fact that Budyko's 
results are yearly mubintoniacten| values, whereas the F)0*% 
values of this study, involving data from 16 October 1973, 
occurs near the time of the oceanic maximum of temperature 
and higher than annual mean cloudiness. This would reduce 
the net flux Fi0* computed here. A latitudinal maximum net 
flux, Fio*, was observed around 40°N both for this study and 
the comparison figures of Malkus. At this location, Figs. 
7(b), 9(b), 10(b) show that there was a minimum opaque cloud 
cover for the time of this study. 

The comparison of outgoing long-wave fluxes, F2 (Table 
II(b)), appears very consistent when comparing only northern 
hemispheric latitudes, considering that Raschke's values are 
based upon the mean for a month rather than that for a single 
day. Here maximum values of very nearly the same order oc- 
cur in the subtropics where minimum cloudiness occurs to im- 
pede the emission of outgoing terrestrial flux from the earth 
and lower atmosphere. There is in both this study and in 
Malkus', a well defined poleward gradient of F2 by latitudes. 

The comparison of the computations of Fio* and F2 made 
Hewervith climatology (Table I1la,b) cannot confirm or deny 
the specific cloud-parameterization formulas (Eq. 2-13a,b) 
used in this study, but they indicate that they have a gen- 


eral reliability. 





(a) (b) 
LAT Fi o* Fi o* Fo F2 


(Malkus) (This Study) (This Study) (Raschke et al.) 


207s . 100 .068 LoSG .40 
10°S 7090 .054 ao 1 439 
0° . 080 .068 OOD 36 
10°N . 090 od .378 .34 
20°N .114 .078 .389 .40 
30 -N 7LO .081 . 388 ou 
40°N 720 LOW .376 34 
50°N lO . 086 . 346 ae 
60°N - 120 . 067 -30D 229 
Global | 
Mean .1018 OL Se . 3601 soouo 
faolewtica). Comparison of net flux at the surface, Fjo*, 


as found by this study for 16 October 1973 and the annual 
mean F,)* derivable from Malkus (1962, Figs. 9,10). 


(b). Comparison of net flux to space, F2, as found 
by this study and by Raschke et al. (1973), based upon 
NIMBUS III heat budget studies. 
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IV. SOLAR RADIATION 


A. COMPOSITION OF SOLAR INSOLATION 

At the top of the atmosphere (k=0) this study assumed a 
Polameconstant of 2.00 ly min7~’ (Joseph, 1971): furthermore, 
this constant was assumed subject to a four percent attenua- 
tion above the tropopause due to ozone and oxygen. Thus the 
effective solar constant at level k=2 in this study is 1.92 
inyemin-* . 

To compute the effective solar insolation at the tropo- 


pause the following formula was used 


m@2e= Sf} - Cos Z (4-1) 
m 
where S = effective solar constant at level k=2 
Cos Z = cosine of the zenith angle, a function of 
the Julian date 
mr = ratio of the actual earth-sun distance to 


the mean earth-sun distance, a function of 
the Julian date. 
The Smithsonian Meteorological Tables (List, 1958) list the 
age b 1O r/r ion howOcCcthober at OOOOGNT as 0.997165; This varue 
was utilized for all four data lines. The cosine of the 


zenith angle was determined by 


Cos Z = Sin » Sin 6 + Cos ¢ Cos 6 Cos h (4-2) 
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where 


@¢ = latitude 

pea=ecolar declination for 15 October = =8° 13.5% 
(Table 169, List, 1958) 

h = hour angle 


The value of sine » was calculated using one of two differ- 
ent formulae, depending on the data-line used for the compu- 
wrerons, in terms of the FNWC map coordinates (I,J) as in 
Eq. (4-3a,b), conversely for these lines one may solve for I 


in terms of Sine ¢ as in (4-3c,d): 


. . en 262 22 See) 2 
Lines 1,3,4 Sin @ 973.752 + 3(32-1)2 (4-3a) 
ee. Omen Oo Coe) . 
ILysLoKey Sin ¢ 073.752 + (32-1)2 (4-3b) 
: _ Cos _ 
Lines 1,3,4 I = 32 - 22.065 [FZ ra g! (4-3c) 
; 7 Cos 6 a 
Line 2 i= S32e— oh. 209 eeenrcar) (4-3d) 
es a RCTS for Lines 1,2 
J = une meee for Line 3 
Go,. <4 50 for Line 4 


Hememin is the abscissa distance on the FNWC grid, Fig. 2, 
and varies by line as described in Sec. II. The soundings 
for lines 1, 2, and 3 were all taken at OOOOGMT with the so- 
lar noon existing at the 180° meridian, therefore the hour 
mimivesetor these three lines were 55°, 10°, and 35°, re- 
spectively. For line 4, the soundings were taken 12 hours 
earlier with solar noon at the Greenwich meridian, giving an 


Houreanele for line 4 of 35°. 
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The composition of solar insolation was partitioned in 
ibs suudy, in accordance with Joseph (1971), by dividing 
Piominseolataon, F(2) at level 2, into two parts, ome of 
which was subject to water vapor absorption but not to Ray- 
leigh scattering, F(A), and one which was only subject to 
Rayleigh scattering, F(S), in a clear moist atmosphere. 


These partitions may be easily computed as follows: 


F(A) 0.349 F(2) (4-4) 


FCS) OFGot FC) (4-5) 


Cloud-reflectivity is inherent in both F(A) and F(S) solar 
energy portions with the presence of two cloud decks. 

To further define the solar insolation partitions, the 
F(A) energy consists of those waveiengths X > 0.9 um where 
absorption by water vapor and carbon dioxide bands is the 
dominant attenuation process for clear air. The F(S) energy 
is composed of all wavelengths » < 0.9 um where water vapor 


absorption is considered negligible. 


B. DISPOSITION OF F(S) INSOLATION 

As described by Jenks (1974), there are two reflective- 
type parameters which must be considered when studying the 
energy contained in the shorter wavelengths, F(S). The 
first parameter is the Rayleigh clear-sky albedo, (after 
Coulson, 1959), which was empirically formulated using the 


least squares best-fit technique by Joseph (1971) to be 
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Mmeaye— .0sa + ~~ 20074 (log, » sec Z)] 


ee 
=“ (4-6) 


Boe wl OLo. 20 mb 
Since mean sea level is very nearly 1000 mb, this study as- 
sumed that 1/Pp9~1. The value of Sec. Z is the inverse of 
Eq. (4-2). 

The second reflective neeeniC ces considered is the surface 
albedo, which is a function of the cosine at the zenith 
angle over ocean surfaces. The surface albedo, as utilized 


by this study, follows Gates et al. (1971). 
mcGyme— max 1.06, .06 + .54 (.7 = Cos Z)} (4-7) 


As described in Sec. III, four distinct cloud cases are pos- 
Sible with a two-layer cloud model. The disposition of F(S) 
under each of these cases will be discussed in the remainder 
of this subsection. 
imewecrear Sky Case 

With clear sky the F(S) insolation is modified only 
by Rayleigh sky reflectivity, a(R), and surface reflectivity, 
a(G), as given by Eqs. (4-6) and (4-7). nenice (1974) points 
out that with multiple reflections possible between the 
earth and atmosphere, the clear sky insolation absorbed at 


the earth's surface (k=10) after scattering is given by 
Gia e ) = F(S)[1-o(R))[t-o(G))7[i-a(R)a(G)] (4-8) 


where the notation (0,0) denotes CL(1) = CL(2) = 0.0 as dis- 


cussed in section III. 


43 





2. Cloudy Sky Cases 


With the presence of cloud layers, an additional in- 
So0lation modifier must be considered, cloud reflectivity. 

In this study the values suggested by C.D. Rodgers (1967) 
were adopted, namely, R(1) = 0.54 for the higher clouds (be- 
tween k=4 and k=6) and R(2) = 0.66 for the lower clouds (be- 
tween k=8 and k=9). 

Consider first the case (1,0) where the fraction of 
the grid area has overcast clouds only at the high level, 
that is where CL(1) = 1.0, CL(2) = 0.0, and R(2) = 0.0 So 
that the earth-absorbed F(S) insolation is 

ee GeO) Cs)Ci-h(1))C1l—-a(G))/(1-R(1)a(G)) Cao) 

The F(S) insolation absorbed by the surface above 
which there iS present only a lower cloud overcast, [denoted 
PwaecOe) | where CL(1) = 0.0, CL(2) = 1.0, and R(1) = 0.0 is 
given by 

1S10(0,1) = F(S)(1-R(2) ) (1-0 (G) )/(1-R(2)a(G) ) Cao) 


If overcast clouds are present in both layers, then 
the conditions CL(1) = CL(2) = 1.0 [or in abbreviated symbol- 
ism, case (1,1)] applied, and reflectance also exists between 


the cloud layers as well as between clouds and space. Thus 


the absorbed FCS) insolation at the earth's surface @s for- 


mipvated aiter Arakawa (1972) 1S given by 
PomlO@in dt )je= F(S)[(1-R(1))CI-R(2) (1-0 (G)))/(1-CCRC I) RO) 


+R(2)a(G)+R(1)a(G) )+2(R(1)R(2)a(G) )] 45) 
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a formula which resembles very closely a combination of Eqs. 
(4-9) and (4-10). 

Thus far in this subsection the sky was considered 
either overcast or clear at each level without regard for 
meeactual cloud conditions CL(1) and CL(2). 

3. Composite F(S) Insolation 

If the areal coverage weights, as defined by Eq. 
(3-4), are applied to the appropriate right sides of Eqs. 
(4-8), (4-9), (4-10), and (4-11), then the composite surface- 
absorbed insolation, IS10, results 


IS10 = 1510(0,0)W(0,0)+1S10(1,0)W(1,0) 


| (4-12) 
Sebo toco, t)wCO,1)+1S10(1,1)W(1,1) 


It 1S apparent that F(S) insolation reflected to space may 


be found by subtracting IS10 from F(S). 


C. DISPOSITION OF F(A) INSOLATION 

This subsection deals with that portion of the solar in- 
solation which is absorbed by atmospheric water vapor and 
carbon dioxide, namely the longer wavelengths. The discus- 
Sion will again be divided into the four different cloud- 
layer cases. 


1. The Clear Sky Case 





Formulation enabling the determination of absorption 
values for particular layers of interest was provided by use 
Steiae Manabe-Moller absorptivity function which may be gen- 
eralized as follows 


meee = 9. 271(U(2,2k) Sec 2]: °°? CARES) 
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Hese absorptivity a is applied over the pressure-scaled water 
Mapeomemass between levels 2 and 2k (Fig. 1) along the slant 
Poavneelt the zenith angle Z. Recalling that the water vapor 
mass above level 2 was assumed negligible, the absorbed solar 


misetaction, tor example, in the layer (2,6) is written as 
Mee oO more 7 FCA)(UC2,6) Sec Z)° 7°? (4-14) 


and the absorbed solar insolation in the layer (6,10) is 


given by 
Neo tO) = sA(2,10) - A(2,6) (4-15) 


For a schematic depiction of the F(A) insolation disposition 
in the clear sky case [see Fig. 3 in Jenks (1974)]. 

Homepuaita that portion of FCA) insolation which is 
absorbed at the earth's surface, the atmospheric transmissiv- 
mvt acc,10)) and the surface transmissivity after reflec- 
tance are combined (by the product-law of transmissivities) 


to give IA10 as follows: 


epee — 9h (A){1—-.271(/U(2,10) Sec Z)-?°*}(1-a(G)) 
(4-16) 


Here a(G) is defined by Eq. (4-7) 

In the remainder of this subsection, which discusses 
the cloudy layer cases, representative cloud reflectivities 
‘and iene absorptivities wWeremcopred, alter C.D, Rodgers 
“Moo werent the two possible cloud layers. The cloud re- 
flectivities were RA(1) = .46 and RA(2) = .50 while the 
cloud absorptivities were A(1) = .20 and A(2) = .30. Note 


Piet eeme tet lectivities for the F(A) wavelengths differ from 
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those adopted for the F(S) wavelengths. The procedure of 
considering the cloud conditions to be overcast whenever 
they appear and then applying the appropriate weighting fac- 
tors in the composite summation will again be followed as in 
mec. IV.B. 

To simplify the discussion for the elend=couenna 
cases, the (1,1) case will be presented first as it contains 
representative type equations for the remaining two cases, 
(1,0) and (0,1). 

2. Overcast in Both High and Low Cloud Layers 

The following set of formulas detail the disposition 

of incoming solar insolation (F(A)) from level k=2 to the 


earth's surface to permit determination of the amount of in- 


solation absorbed by the atmospheric layers and by the earth's 


surface. The dashed Separation lines are introduced to sub- 
divide the absorption and reflection physics of the model 
into subsections which permit the analysis to proceed more 


or less within successive 200 mb layers (see also Fig. 3). 


feel 7 1 | Uieee4 ) Sec Z)> 7°) 


F4+ = F4V(RA(C1)) 


(4-17a) 
Pwetie—en(A),271,002,4) Sec Z1-*°° 
Pye 4+. 971(U02,4) Sec Zy°°?°* 
F64 = F44(1-RA(1)-A(1)) 
A46 = F44(A(1)) (4-17b) 


A26 = A24+¥ + A24t+ + A46 
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Figure 3. Schematic representation of F(A) insolation 


disposition in the case of two overcast layers. 
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Mies =—1—=.271[(U(6,8)5/3]- ?°3 


ipslor= 1-271) U9, 10)5/3)" 2° 3 


F8) = F6\(TD68) 
F8t = F8+(RA(2)) 


F6t = F8+(TD68) 


F6)4 = F6+(RA(1)) 
F844 = F644 (TD68) 
AGS) = FG) - F84 
AG8t = F8t - FG6t 


A684) = FGYY - F8+4 


A68 = A68Y + A68t + AB8YY 


F104 = F9+4(TD910) 


ry 
}4 
© 
BS 
ll 


F104 (a (G) ) 


F9+ = F10t(TD910) 


ry 
© 
Ea 
ae 
lI 


F9+ (RA(2) ) 

F104) = F9¥+¥(TD9I1O) 
A910 = F9Y - FIO 
A910t oop i Ot 


AQ10++ = F9OVY - FILOVY 


A910 Romo, + AGLO+ + AOMOY ¥ 


A810 = A89 + A910 


a a a 
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(4-18) 


(4-19a) 


(4-19c) 








Mmehespreeeaing equations the symbol “A” followea@ 
byea Set of digits refers to the F(A) insolation absorbed in 
the indicated layer, for example, A26 indicates the absorp- 
tion in the layer k=2 to k=6. The symbol "F" followed by a 
Pevallue indicates the amount of insolation passing through 
the particular k level. The vertical arrow modifiers indi- 
eave the direction of insolation passage, i.e., + downward 
insolation, + upward reflected insolation, and ++ downward 
reflected insolation. Terms involving the symbols "TD" as 
expressed by the functions of Eq. (4-18), indicate the 
Manabe-Moller transmissivities for insolation beneath a 
cloud layer. In the latter situation, the term Sec Z@ in 
Eq. (4-13) is effectively replaced by the mean slant-path 
Sec Z = 5/3 (Katayama, 1966). 

Note that the effect of multiple reflections between 
clouds or between the earth's surface and the lower cloud 
has been incorporated to the extent of two reflections, with 
the lowermost reflecting surface absorbing the remaining im- 
pinging insolation. Computations indicated that the insola- 
tion remaining after two reflections was too small to warrant 
further reflections. Also insolation reflected upward from 
a lower interface (cloud or ground) to the base of an upper 
cGloud deck has not been subjected to cloud absorption. This 
tends to reduce very slightly the secondary cloud-absorption. 

Note finally that between any two reflecting surfaces 
there are three contributions to absorption of solar insola- 


Pwonewithin the layer. These contributions are denoted by 
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the arrow-identifiers which indicate the segment of path 
being traversed. 
From Eqs. (4-19c) the impinging F(A) insolation at 


the earth's surface can be expressed 
TRANA(1,1) = F104 + F10+4¥ (4-20) 


The F(A) insolation which is actually absorbed by the earth's 


surface (see Fig. 3) may be written as 
TA10(1,1) = F104(1-a0(G)) + F1OVY (4-21) 


3. Disposition of F(A) Insolation with an Upper Overcast 


With a cloud layer present only in the upper layer 
the equations depicting the disposition of incoming insola- 


tion becomes a Simplified subset of the previous case: 


BGA) S271 U(2,4) See 2° ?*? 


A244 


A24t Penile. 2711U(2,4) Sec Z]°**?} RAC) (4-220) 
eeonilicoyd ) SecwZ)| °°: } 


A24 = A24} + A24t 


wee wee mm mm i i Se i TF ll 


Hete—=T( 8 11-.27110U(2,4) Sec Z)* °°? }11-RA(1)=AQ)) 


> 

— 

oO) 
II 


Hew i= .2711U(2.4) See Zi >?" TAC) (4-22b) 


A26 = A24 + A46 


— jo o9 ooo oom oom oom oo gt wwe iwi we we we ee CTT CTT CC LTCC CC Tae 
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F10¥ = F64(TD610) 
F10+ = F104(a(G)) 
F6t+ = F10+(TD610) 


Hee = FOt(RAC1) ) 


ELOY) = F6+44+(TD610) (4-22c) 
Belov = FEY — F104 
A610+ = F10+ - F6t 


A610++ = F6+¥ - F1OVY 


A610 = A610+ + A610¥ + A6G1O+¥ 


The variables used above are defined in a similar manner to 
mo@sce sin the (1,1) case. The impinging insolation at the 


earth's surface may be formulated as follows 
TRANA(1,0) = F10Y + FI10O+V¥ (4-23) 


while the F(A) insolation absorbed by the earth in this case 


(Fig. 4) is given by 


IA10(1,0) = F104(1-a(G)) + F1044 (4-24) 


foeoesposttion of F(A) Insolation with a Low Overcast 


With an overcast lower cloud layer the disposition 


fone seare as follows: 
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F(A) 

2 F2T 

A24t A24% 
4 F4¢ RA(1 | 

a ee ee 
RA(1) 
FOLLY Fot F64 
A610144 A610 t A610 4 
8 
ie F108 a(G) 
oma fa, 7 ff ff fF Y ff 4 7+ 
1A10 
Figure 4. Schematic representation at F(A) insolation dis- 


position with an upper overcast layer only. 





Maore=shCA).271(U(2,6) Sec Z)-*°? 

eee err ii| UC2,8) See Z)"°**-(_UC2,6) See Zz) ~~ ?} 
Bae HCA) ( i=. 271, UC2,8) Sec Z): 3%?) 

F8t = F8y[RA(2)] 

F6t = F8t{1-.271[U(6,8) Sec Z]*?°3} (4-254) 
F2+ = F8t{1-.271[U(2,8) Sec Z]° 2°} 

A68t = F8t - F6t 

A26+ = F6t - F2t 


A26 = A26y + A26t 


FOY = F8¥[1-RA(2) - A(2)] 


(4-25b) 
A89 = F8+4[A(2)] 
F104 = F9+4(TD910) 
F10t = F10+(a(G)) 
F9+ = F10+(TD910) 
FOV) = F9t(RA(2)) 
F1I0++ = F9¥4(TD910) (4-25) 


AQI0O+ = F9Y - FIO} 
Soot = F1Ot — F9t 
AQ10++) = FOV) - FIOVY 


A910 


Rolo, + AQLOt- + AQLOY + 


A610 A68+ + A68t + A889 + AYIO 


Neainethe variables are defined as stated before in the 
(1,1) case. Likewise the incident flux is defined at the 


Satemos) SUrface as 


TRANA(O,1) = F10+¥ + F1IOV+¥ (4-26) 
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while that portion which is absorbed by the earth, (Fig. 5) 
is 


TALTO(O,1) = FLOV(1-a0(G)) + F1O+4 (4-27) 


9. Weighted F(A) Layer-Absorptions and Surface-Absorp- 


tion Insolation 

As described in Section IV.B.3 on the composite F(S) 
insolation, the values of the absorbed F(A) insolation in 
the layers and of the earth's surface which were calculated 
hOneLhe four cloud cases (0,0), (1,0), (0,1), and (1,1) must 
be appropriately weighted by the values from Eq. (3-4) to 


obtain the composite F(A) layer- and surface-absorbed inso- 


lation. 
A26 = Gm oo eo od oy 12 2 =e (oy 
(4-28) 
oy iy) S23) 
A610 = AG10(4 9 yW#(0,0)FAG10/4 gyW(1,0)FAG10/ 4 4 NCO, 1) 
, (4-29) 
FA610(4 4)¥C4,1) 


IA10 = IA10(0,0)W(0,0)+IA10(1, 0)W(1,0)+IA10(0,1)W(0,1) 
(4-30) 
+JA10@@,1)w(1,1) 


6. <Absorptivity (ABA) by Layers 


It was convenient to work with the absorptivity as 
well as the actual insolation energies absorbed by each lay- 
er and at the surface. To obtain the absorptivity, the 
total incoming insolation at k=0 must be known. This value 


called FADJ is formulated as follows 
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Figure 5. Schematic representation of F(A) insolation 


disposition with a lower overcast layer only. 


OG 





FADJ = 2.00 cin)” Cos Z (4-31) 


By taking the ratio of the insolation absorbed in the tropo- 
Spnere to that incident at the top of the atmosphere, the 
Mir tropospheric absorptivity is formed 


A26 + A610 


ABA = FADJ | (4-32) 


Similar formulas may be expressed for the partial absorptiv- 


eres for the layers (2,6) and (6,10). 


D. REFLECTION OF F(A) AND F(S) ENERGY 
To determine the albedo and fractional absorptivity of 
the earth's surface, both the F(A) and F(S) fractions of so- 
lar insolation must be considered. 
1. Albedo 
The amount of F(A) insolation which is reflected 


back to space is determined by the following 

REFA = F(A) - A26 - A610 - IA10 . (4-33) 
here A26, AG10, and IA10O are defined by Eqs. (4-28), (4-29), 
and (4-30), respectively. 

Because the F(S) wavelengths do not undergo atmos- 
pheric absorption, the reflected F(S) insolation was easier 
eo Calculate 


Reps = FCS) -— 1810 (4-34) 


IS10 was defined by Eq. (4-12). By combining Eqs. (4-33) 
and (4-34) the total reflected insolation to space is de- 


fined ‘by 
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REF = REFA + REFS (4-35) 


The resulting REF values were calculated at each gridpoint 
and are meridionally displayed by latitudinal interpolation 
mmprigs. V(a),...,l0(b). 

Bye waking the ratio of the total refbectedvenerzy to 
the total incoming energy the earth-atmospheric system albedo 


is defined 

ALB = REF/FADJ (4-36) 
where FADJ is defined by Eq. (4-31). 

2. Composite Absorptivity (ABG) by the Earth's Surface 

The fractional absorptivity of the earth's surface 
was determined by the ratio of the total insolation absorbed 
by the earth to the total incoming insolation. The total 
insolation absorbed by the surface is determined by the prop- 
erly weighted sums of IA10 and IS10, defined by Eqs. (4-30) 


and (4-12), respectively. 
BEGe— ClIA1O + IS510)/FADJ (4-37) 


The result of Eq. (4-37) gives the weighted absorptivity of 
earth. It is also of interest to determine the total inso- 
lation incident at the earth's surface just before absorp- 
tion by earth. The result is given by the properly weighted 


SlummOoneberns OL the form 


TRAN = TRANA + [IS10/(1-a(G))] (4-38) 





Pieshould be noted that each of the four terms of form TFSi0 
mea. 4-12) had a factor (1-a(G)) in the numerator so that the 
transmitted F(S) insolation incident at earth before absorp- 
tion need only be divided by (1-a(G)). The parameter TRAN 

is a simulation of the insolational energy incident at a 


pyrheliometer. Note finally that the ratio 
ATRAN = TRAN/FADJ (4-39) 


represents the transmissivity of the troposphere. 
3. Computational Check 
Recalling the assumption of a four percent attenua- 
tion of the incoming solar insolation in the stratosphere, 
the summation of the fractional absorptivities ALB, ABA, and 
ABG should equal 0.96. This value provides a good computa- 
tional check for this model at each gridpoint and was in 


fact satisfied over the entire grid. 


E. STATISTICAL TESTS 

To further verify the foregoing calculations, statisti- 
cal computations were performed on several key parameters 
utilizing linear regression techniques from the BMD set of 
Pruiatteale computer programs (Dixon, 1973). 

1. Clear Sky Case Regressions 

Three linear regressions were computed for the clear 

sky case, the first of which was the clear sky albedo, 
ALB(0O,0,Z2). The best-fit equation, using the secant of the 


genith angle, Sec Z, as the predictor, was computed to be 


ABC O,0,Z) = 0.00650 + 0.09052 Sec 2 (4-40) 
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With a regression coefficient R = .964. This highly «sagni f- 
teame CcOrrelation coefficient confirms the model results 
witch imdicate that Sec Z is a good predictor of clear-sky 
wpeao®. in this connection, it should be recalled that a(G) 
and q«(R) were modeled to be functionally related to Sec Z. 


The mean global clear-sky albedo was 
ALB(O,0) = 0.164 


which corresponds to the mean secant of the zenith angle, 
Sec Z = 1.743, for this study. The statistics show that the 
larsest contributions to error in the fit by Eq. (4-40) oc- 
curred at largest zenith angles, which occur at high lati- 
tudes. 

The two other clear-sky linear regressions performed 
were ABA(O,0,M) and ATRAN(O,0,M), where the predictor, M, 


was the water-vapor mass path length given by 
a. a 
fe—eaiecee Z)-( log,,(U Sec Z)-] (4-41) 


This parameterization of the water-vapor mass is formulated 
Meter hanson (1967, 1971), who used a similar sheer beusmee 
predictor in his empirical equations for ABA with both clear 
and partly cloudy skies. The resulting clear-sky regression 


equations in this study for ABA and ATRAN were 
ABA(O,0,M) = 0.10320 + 0.05530M (4-42) 
with R = .9925 and 


ATRAN(0,0,M) = 0.77494 - 0.04823M (4-43) 
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with R = .7142. The sample means of the statistics which 


entered into (4-42) and (4-43) are 
ABA(0O,0,M) = 0.146 


ATRAN(0,0,M) = 0.738 


ee a. 
M = 0.776 (gm cm-? )? 


based on the original 93 cases. 
2. Cloudy-Sky Case Regressions 

In addition to the clear-sky fractional amounts of 
ALB, ABA, and ATRAN, linear regressions were formulated for 
the composite cloudy-sky case, i.e. ALB[CL(1),CL(2)], 
ABA[CL(1),CL(2)], and ATRAN[CL(1),CL(2)], respectively. The 
regression relationships between the clear- and cloudy-sky 
cases were determined for the same zenith angle at gridpoint 
mero oewsing the total opaque cloud cover (CL) as the predic- 
tor. The total opaque cloud cover approximates the effective 
cloud cover existing in a two-layer cloud model and iS ex- 


pressed as 
Glee ck( lj + CL(2) -— CLCL)CL(2) (4-44) 


The relationship between the cloudy-sky albedo and 
the corresponding clear-sky calculation was obtained by 


least squares with the following results 


ALB[CL(1),CL(2)] = ALB(0,0)[1+2.09929CL+0. 84608CL? ] (4345) 
and R = .9292. The means of the parameters used in Eq. (4-45) 
aie 
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MieteLOM. chi2)) = 07406 


ALB(0O,0) = 0.164 


Clie 107662 


The dependence of Eq. (4-45) upon CL is indicated by 
the high correlation coefficient (R=.9292). The global mean 
albedo for this study of 0.406 is somewhat higher than that 
reported in recent literature (ALB=.300) by Raschke et al. 
(1973). However, Raschke's figure is an annual mean which 
includes the somewhat smaller cloud amounts over land. 

The tropospheric absorptivity of air in the composite 
cloudy case was developed analogously to (4-45) by least 


squares and gave the form 
MBemCi), CL(2)] = ABACO,0)[1+40. 22305CL+0. 17261CL*] (4-46) 
with R = .927 and mean statistics for the above parameters 
ABA[CL(1),CL(2)] = 0.180 
ABA(0,0) = 0.146 
CL = 0.662 


The direct relationship between solar absorptivity and cloud 
cover is in general agreement with the results of Plante 
(1972) who found increased absorptivity with increasing CL. 
This relationship could be the reason for the recent reduced 
eoimabes ol global albedo [Raschke et al. (1973), Vonder 
Haar and Suomi (1971)]. The significant value of R indi- 
cates the degree of sensitivity between cloud amounts and 


atmospheric absorptivity. 
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The final statistical regression equation formulated 
deals with the cloud-sky transmissivity of the troposphere 
Betative LO the clear=sky calculation. A@ain te tote 


Opaque cloud cover has been used as the independent variable 


ATRAN[CL(1),CL(2)] = ATRAN(0,0){1-0.51477CL-0.14646CL? ] 
(4-47) 


Here the correlation coefficient R = .997 and the mean sta- 


tistics over the 93 sample soundings are 
ATRAN{CL(1),CL(2)] = 0.543 
ATRAN(0,0) = 0.738 
CL = 0.662 


Note the high degree of significance, (R = .997) of Eq. 
(4-47), which formulation agrees with empirical results of 
Quinn (1971) and numerous other workers cited by Quinn in 
his 1971 paper. 

In general, the solar-model statistical relation- 
Ships developed in this study show consistency with empiri- 
cally reported observations, which have already been noted, 
when consideration is taken for seasonal effects and the 


latitudinal limits of this study. 





V. SENSIBLE AND LATENT HEAT TRANSPORT 


AT THE SEA-AIR INTERFACE 


A. GENERAL CONSIDERATIONS 

For a complete discussion of the heat balances of the 
earth and atmosphere for 16 October 1973, the vertical 
transports of sensible and of latent heats across the sea- 
air interface must be considered. For this purpose, the 
turbulent flux models for the planetary boundary layer, as 
used in the FNWC primitive equation model as described by 
Kesel and Winninghoff (1972) and Kaitala (1974) , was adapted 


here virtually without change. 


B. EVAPORATION 
The determination of the surface latent heat flux, E, 
over the ocean is based upon its large scale parameteriza- 


Peon prinecipally after Langlois and@Kwok (1969) 


q_7q 
x 
E = Leyes (A (5-1) 
3. CUS 9-410 
* 
Bie (Feb 


Here L is the latent heat of vaporization given by 
L = 596.73 - .601 Tio (cal/gm) 


inesether parameters in (5-1) are formulated as follows: 


Pio. Pio/Rytio 





00) 2/7 


85 = T5990 


8 = Tio 
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R 9 10 
nee oe 1000 
eee 6 1M “S00 
ex 
4x74, = ,62197 1000 -q( 900) 


Note that at level k=10, the assumption P,, = 1000 mb is con- 
Sidered valid and e, 1s the vapor pressure at the top of the 
constant flux layer. Also q(900) is the mixing ratio at k=9 
as defined by Eq. (2-4). In Eq. (5-1) K* and a*® are turbulent 
transfer coefficients having values 2 x 10° cm* sec™! and 

5 x 10* cm/°K, respectively, according to Kaitala (1974). 
However, values of the eddy coefficient, K*, of 2 x 10° and 

1 x 10° were both tested in this study and it was decided 
primarily to use only the latter value. Furthermore, it 18s 
felt by this author that upon further testing, with a broader 
data base, it may be feasible to tune K* between the two 
values tested so as to make it dependent either upon latitude 


or season, or other pertinent variables. 


C. SENSIBLE HEAT TRANSPORT 

Sensible heating is computed as a function of the tem- 
perature difference between the surface, T;,, and the air 
near the surface, TY which is at thestop of the constant 
flux layer. Since wind values were not carried in this 
study, the technique used by Mintz and Arakawa, as described 
by Kaitala (1974) and by Kesel and Winninghoff (1972), was 


MretratctLempted for sensible heating 
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Kx @.=0., 





eo. Cc (—<—<—— = oe an 
r 10 pe ; = Des =, (5-2) 
+a 
ET 
Here 5 ~eeecouvecal/(pm°K), and, since the constant Miux lay= 


er was assumed to be quite thin, Oa 0s is a eritwreal 
lapse rate (3 x 10°° °K/cm), and the remainder of the terms 
are defined as in the discussion following Eq. (5-1). The 


micrometeorological equivalent of Eq. (5-2) is given by 


Kaitala (1974) as 
Hp — pS sa) x (5-3) 


It can be seen from Eqs. (5-2), (5-3) that the temperature 
lapse, ig for .the thin layer between the surface and the 
top of the constant flux layer becomes very critical in ee 
termining the amount of sensible heating (cooling) for each 
sounding location. 

In Eq. (5-2), the expression within the brace determined 
the sign of the sensible heat flux in the vertical. In ap- 
proximately one out of three cases, the braced expression 
pave rise to negative values of He, ihverneae. the correspond- 
ing latent heat flux, E, of Eq. (5-1) was systematically 
positive. Based on Fig.15 in Malkus (1962), which shows H,.?0, 
imeremscems to be no valid basis for having an erratically 
downward distribution of sensible heat over the open oceans. 
Consequently, this study assumed that E was properly as- 
Seccedeby Eq. (5-1), and that a réasonable set of values of 


Hy. could be obtained by means of 
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H. = rk (5-4) 


Here the Bowen ratio r was determined from a set of best es- 
famaves aS a function of latitude. 

Using 20 values of r and their functionally related lati- 
tudes, taken from published results of Sellers (1965), 
Budyko (1956), and Jacobs (1951), a linear régression was 


best-fitted to [Sin 6] to give 
r = 0.10910 - .41996 |Sin 9] + .89645 |Sin |? (5-5) 


with a correlation coefficient of R = .982. Thus, using 
Eqs. (5-4) and (5-5), the sensible heat transport was ob- 
tained at each gridpoint, while always maintaining proper 


sign sense with E which was virtually always positive. 


D. COMBINED TURBULENT HEAT FLUX ACROSS THE SEA-AIR INTER- 
FACE 
To get a better feel for the effect that sensible heat- 
ing has at the surface, the summed terms E and Hp were uti- 
lized as heat-loss quantities in the total budget at the 
puirace s therefore in terms of large scale “convective dis- 


position 


E + Hp = E(1t+r) (5-6) 


where r ranges in magnitude from approximately .05 to .95 


with increasing latitude. 
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VI. MERIDIONAL CROSS-SECTIONAL DEPICTION OF 


THE HEATING MODEL COMPUTATIONS 


A. GENERAL DESIGN OF HEAT-MODEL OUTPUT 

In this section, the computational subroutines already 
developed (Secs. III, IV, and V) are combined to perform a 
Single time-step of the heating program developed for adapta- 
tion to the FNWC prediction model. The results represent 
computations of the relevant heating functions at each grid- 
pommre ol the four meridians (Fig. 2). 

For the three Pacific cross-sections, the FNWC gridpoint- 
processed analyses for OOOOGMT, 16 October 1973, have been 
used, whereas for the 35°W (Atlantic) meridian, the FNWC 
analyses from 1200GMT, 15 October 1973, were utilized. For 
Simplicity the heating model computations will be referred 
fwomas EhOSe pertaining to 16 October 1973. 

Fig. 6 depicts the key to the computations made at (I,J) 
gridpoint soundings, shown in the form of Table I(b), and in 
effect represents the flow diagram of computations proceed- 
igi romecop—-to-bottom at each gri@point. Figs. 7, 8, 9, 
and 10 show that the gridpoints considered are at 5° lati- 
tude increments, but this has been done merely for simplifi- 
cation of the interpretations and for comparisons with 
climatological data. In order to deduce a heating parameter, 
Smijyeeat a non=integral grid-value J corresponding to lati- 
tude ¢? (in whole multiples of 5°) Sambaer tiga enecub ,vemiaier— 


polation scheme was utilized of the form 
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per) = 2u6t=2)(1-3)(1-4) ,°@, (1-1) (1-3) (1-4) 


(1-2) (1-3)(1-4) (2-1)(2-3)(2-4) 
(6-1) 

peeled )Cl=2)Cied) | 8, Cl-1) G23) 

(3-1)(3-2)(3-4) (4-1) (4-2)(4-3) 


Here I is determined by an equation of the form (4-3c), the 
points 1,2,3, and 4 are any four successive integral values 
of I with I lying in general between points 2 and 3, although 
the subroutine also works well if I is non-centered in the 


I-array. 


B. CROSS-SECTIONAL DATA AT LEVEL k=2 

The terms explained here are those occurring at the lev- 
el k=2. The discussion of all insolational parameters dis- 
cussed previously in section IV pertained to the specific 
time of day corresponding to the hour angle h at the instan- 
taneous times under consideration. Thus the incident inso- 


lation was given by Eq. (4-1), i.e. 


r 


F(2) = S CG) - Coc 
Im 


with Cos Z given by Eq. (4-2). For the specific map times 
Miemlocal solar times gave h=35°, 10°, 55°, amd 35° for 
eress-sections 1,2,3, and 4 depicted in Figs. 7, 8, 9, and 
10 respectively. 

For each specific zenith angle Z corresponding to lati- 
tude ¢6(1,J,t) the absorbed insolations A(276,t)y >and #(6,10,8 ) 
of Eqs. (4-28) and (4-29) have been determined for the given 
time t of 16 October 1973 for which the specific hour angles 
were applicable. The instantaneous parameter Peect) o1 Eq. 


(4-35) was determined by the equation 


Og 








REF(t) = F(2) - A26 - A610 —- (IA10O + IS10) (6-2) 


where the expression in parentheses is identical to the in- 
stantaneous absorbed insolation of earth given by the numer- 
ator of (4-37). 

fhe net terrestrial flux F,* was derived in Eq. (3-11) 
and was not subject to any systematic diurnal variation. 
Thus F,* may be taken as constant for a 24-hour period. 

Note that the first symbol at k=2 in Fig. 6 is not F(2), 
but QAVE, defined for each gridpoint by the analog of (4-1), 
namely 


QAVE = S(—)7?2 Cos Z (6-3) 
Um 








Here Ccs Z is the 24-hour average of Cos Z, and iS given Dy 





Cose4 — 9 (h sind Sin 6 + Cos 6 Cos 6 Sin H]/t (6-4) 
where 
H = are Cos [-Tan 6 Tan 6] (coop 


is the hour angle from local noon to sunset at the gridpoint. 
Thus the instantaneous values of incident insolation and of 
reflected insolation have been converted from F(2) and REF 


to their 24-hour averages denoted by 








QAVE = F(2) Cos Z/Cos Z (6-6) 
and 

QREF = REF Cos Z/Cos Z (6-7) 
respectively. 
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The reason for the conversion in (6-6) and (6-7) is to 
convert the net incoming insolation at k=2 to values that 
are comparable in time-scale to those of F,*. This device 
has the advantage of showing whether the resultant averaged 
heat-transfer quantities are compatible with the expected 
approximation to a heat budget for the given date. However, 
if it were desired to exhibit the expected diurnal variation 
of heat-transfers across the layers (2,6), (6,10) and at the 
ground, the quantity QAVE and all other solar-heating rates 


to be discussed would be multiplied by the factor 
C€os Z/ Cos Z) 


assuming that all the layer absorptivities and the solar re- 
mlectivity are constant during the day for which the diurnal 
variation of heat transfers are of interest. 
With the solar quantities averaged for the day, one ob- 
tains 
BALT = QAVE - (QREF + F,*) (6-8) 


as a 24-hour average radiational balance in cal cm? min7? 


for the tropopause. 


C. CROSS-SECTIONAL RADIATIVE TRANSFERS IN LAYER (2,6) 

At this layer, in Fig. 6, the heat-transfers are assumed, 
as in the FNWC heating model, to be radiative only. For con- 
venient comparison with climatology, the heating (cooling) 


rate 1S given by 


BAL26 = Q26 - F26 (6-9) 
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milere O26 1s the daily average absorption in layer (2,6) and 
is defined relative to A26 by an equation analogous to Eq. 
(6-6). The terrestrial cooling rate F26 is defined by Eq. 
(so-lZa). 

Note that all values of BAL26 were negative in all cross- 


Bectwonal diagrams, Figs. 7,...,i0. 


D. CROSS-SECTIONAL DATA FOR LAYER (6,10) 
1. Radiative Transfers Only | 
Neglecting temporarily the turbulent transfers at 
the earth-air interface, the present heating model computes 
the 24-hour average radiative cooling computed by the differ- 


ence equation 
BAL610 = Q610 - F610 (6-10) 


in a manner analogous to Eq. (6-9). The parameter BAL610 
Voiomiesautyve at all gridpoints of Figs. 7,...,10. In fact 
the ratio R = F610/Q610 varied with latitude as follows 
OG 10° Z0q 30° 40° DOG 60° 

Reo 2 Lome oie 2.695 3.912 4.122 4.903 
Even when the instantaneous (near-noon) radiation balance 
was considered, the difference remained at least slightly 
negative in low latitudes corresponding to cooling. In high- 
er latitudes, instantaneous cooling rates were never reversed 
by consideration of the instantaneous computation of the ra- 


Giationmeablance in Layer (6,10). 


a 





2. Heat Balance at the Earth's Surface 
In Section V.B, the earth-to-air turbulent transfer 
of latent heat, E, has been modeled after Kaitala (1974), 


through Eq. (5-1). The sensible heat transfer, H. given by 


iH, = rk 


where r is the Bowen ratio determined empirically by Eq. 
(5-5) as a function of latitude. 
The sum of the turbulent heat transfers, HL +E, has 


been subtracted from the average radiation at the earth's 


surface 
R = QABG —- Fj, 0* (6-11) 
to form the average warming rate, BALB, at the surface 


BALB = QABG —- Fio9* - (H, +E) (6-12) 


where QABG is related to ABG of Eq. (4-37) by an analog to 

Eq. (6-6) and F,;,9* is defined by Eq. (3-10) which is taken 

as constant for a 24-hour period as was the case with F,*. 
The parameter BALB has been shown.in each cross-sec- 

tion as a function first of gridpoint on each meridian, and 

Subsequently interpolated to 5° latitude increments as indi- 

cated in connection with Eq. (6-1). 

| 3. Sensible Heat Modification of the Atmospheric Columns 

In Figs. 7,...,10 the sensible heat transfer, Hp, 

which crosses the earth-air interface has been added to the 


Piimeon witimeicmmweyer (5,10), (cil. Kaitala, 1974). In this 


connection, the net radiative loss in the layer (6,10) has 
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been subdivided equally into two halves, BAL68 and BAL810 
(excluding H.), equal to half the radiative loss BAL610. 

The sensible heat transfer He, which has been allowed 
Gemcditiuse Onlyeinto layer (S,10), is shown as item (n) of 
ever (8,10) in Fie. 6G. In the numerical computations, H, 
is the smallest of all the heat-transfers Bebe dered ii eany 


aemospheric column presented in Figs. 7,...,10. 


E. MERIDIONAL CROSS-SECTIONS OF THE VERTICAL HEAT BUDGET 

Note that the single-time set of heating computations 
Houmecach @©l Che four meridians, 125°W, 170°W, 145°E, and 
35°W, are shown in Figs. 7,...,10 respectively. These fig- 
ures show the individual computations of the warming and 
cooling which are schematically formatted in Fig. 6, together 
with the layers to which they contribute in the heating mod- 
eo). 

Note also that the four meridional cross-sections de- 
picted each have two sections: section (a) depicts Peon ES 
for a “tropical™ section with gridpoints extending from 20°S 
(or the southernmost latitude to which the meridian in ques- 
Plonme<tendsjto 2o N, while Section (b) of Figs. 7,...,10 
depicts the gridpoint computations (averaged for 24-hours) 
Sverre uhemlaLioudinal range SON to 65°N (or the northernmost 
latitude to which the meridian in question extends). Recall 
the latitudinal boundaries for the selected meridians were 
required to remain over the ice-free ocean on the FNWC grid. 

The latitudinal variation of both the short- and long- 


wave calculations, as observed in Figs. 7,...,10, seem to 
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Haver unttorn @radients across each meridian in spite of the 
differing climatic zones represented. Hence it was reason- 
able to formulate a mean radiational balance model for 16 


October 1973, which is presented in Section VII. 


io 
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125°W Longitudinal cross-sectiion, tropical section. 


Values computed from data for 16 October 1973. 
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Refer to Fig. 


145°E Longitudinal cross-section, tropical section. 


Values computed from data for 16 October 1973. 
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145°E Longitudinal cross-section, higher latitude section. 


Figure 9(b). 


Values computed from data for 16 October 1973. 
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Refer to Fig. 


tropical section. 


35°W Longitudinal cross-section, 


Figure 10(a). 


Values computed from data for 15 October 1973. 
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Refer to 


35°W Longitudinal cross-section, higher latitude section. 


Figure 10(b). 


Values computed from data for 15 October 1973. 
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VII. MEAN BUDGETARY BALANCE MODEL 


By averaging key parameters, discussed in Section VI, 
across the four possible meridional lines at each 5° latitude 
increment from 0° to 65°N, a mean heat budget model was con- 
structed which may be considered representative of Fall in 
the annual cycle. The advantage of such averaging is to fa- 
cilitate the comparisons of this study's results with other 
published heat budget results (cf. Budyko, 1956) even though 
such results are generally annual means. Also, by averaging 
over the four meridians studied here, meridional differences 
tend to be smoothed out, so that the zonally-averaged heat 
budgets of this study may be compared with others in the lit- 


erature. 


A. RADIATION MODEL OF THE OCEAN-ATMOSPHERE SYSTEM 

Fig. 1l1l¢(a) depicts the schematic format of the mean ra- 
Gatlon meriaional cross-section (Figs. 11(b) and (c)). 
Turbulent transfer parameters have been omitted from this 
display and only Northern Hemisphere values are given. 

Note that the radiational quantities for the layers 
(6,8) and (8,10) of Figs. 7,...,10 have been combined at 
temic iMebattuudes om Fig. Ll as for the layer (6,10) and 
that all resultant values are meridionally averaged where 
possible. Since only the 170°W meridian contains results at 
G5°N, the value at this latitude is not representative of the 


zZOnal-averaging process described above. 
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Mowtacilitate a ygrapnrical display of the results con- 
tained in Fig. 11 and for ease of comparison with published 
fo eect ania hogied lereceutts {(Malkus, 1962), the model dis- 
played in Fig. 11 has been used as the basis for testing the 
following component heat balances: (a) the earth-atmosphere 
system balance (using symbols consistent with Malkus (1962)), 
(b) the atmospheric radiative balance, and (c) the earth- 


Surface balance. 


B. RADIATION BALANCES OF THE OCEAN-TROPOSPHERE SYSTEM 
imeem, Pam. 12> consists of a radiation bal- 
ance composed of three parameters, Ro; Ro and R, which are 
defined (after Malkus, 1962) as follows: (1) Ro 1S mean ra- 
Gdiative energy gain aoc) rate of the ocean-troposphere 
system which was referred to as BALT at k=2 in Fig. 11(a); 
(2) Ro is the mean radiative cooling rate in the troposphere 
which is defined as the sum of BAL26 and BAL610 from Fig. 
iia) and (3) R, the mean radiative warming (cooling) rate 
at the earth's surface which was identified as BAL1O in Fig. 


li¢a). The relationship between the three radiative param- 


eters can be expressed as 


Re = kh + Ro ci. 


Table III lists quantitative comparisons of Ro; Ro: and 
Doronila tis cuuady tor 16 October 1973, with values 
reported in Malkus (1962) for the climatological year, and 
ValieommoporlLedsby eenks (1974) for 25 April 1973. The lat- 
itudinally-averaged values are obtained by weighting each 


value by the coSine of the latitude 


oo 
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Figure 12. Radiational balance at the tropopause (Ro); in 
the atmospheric column between k=2 and k=10 (Ro); and at the 
Simercem in) for 16 October 1973. 
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Table Iii. Comparisons of Radiational Balances (a) at the 


tropopause Re: Cen the vertical column Ro between k=2 and 


k=10, and (c) at the earth's surface R of the meridionally 


averaged values as found in this study denoted by subscript 


ue la ’ 


ork 


The corresponding annual values reported by Malkus 
ClIGZyware denoted by subscript "2", 
i oeOunamoyelenks (fo7d) denoted by subscript "3". 


VMeseare in ly min? , 


diceehiose tor 25 April 


All 
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It can be noted from Table III that although the Be val- 
ues found by this study are consistently lower than those re- 
ported by Malkus for the year, the seasonally higher values 
reported by Jenks tend to offset this difference. The same 
holds true for the balance at the surface (R). However, 
both seasonal studies show greater radiative cooling in the 
Pea aNe ee (Ro) than the yearly climatological values which 
evidently was obtained after Budyko (1956) as a residual be- 
tween Ro and R, as functions of latitude. It is noteworthy 
that Jenks' (1974) Re Beswrese(CCGlumn (ec) Of Table Il!) ere 
greater than those obtained here by .094 ly min~! (in the 


mean), and of this surplus the additional heating rate 


De OMmOstelyvemine is found at the surface. 


C. TROPOSPHERIC HEATING 

Fig. 13 graphically depicts the mean tropospheric heat 
budget as found in this study. The parameters presented in- 
clude: Ro which was discussed in Sec. VII.B, (E+H, ) which is 
the turbulent heat transfer rate due to latent and sensible 
heat transport realized in the atmosphere, Oe. which is the 
net heat flux divergence rate in the troposphere required to 
bring about a heat balance (which would not be expected to 
be achieved for the time scale here), and So the storage 


term. These terms are related by 


OZ 
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Table IV. Comparisons of trapospheric radiative cooling rate 
(R,), turbulent transport warming rate (EtH,), and net trop- 
ospheric heat loss rate Oe as found by this study for 
Po OCwoObecr slo /s denovea by subscript "1"' and as reported by 
Malkus (1962) for yearly climatology denoted by subscript "2". 
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The storage rate Sa has not been separated specifically from 


Q 


ya? Dut over the annual cycle it is assumed zero. Assuming, 
however that Oe Oe if the weighted average values from 
Table IV are set into Eq. (7-3), to determine the mean stor- 


age cooling rate, Sa? the latter turns out to be 


I 


Se = 1 


- eae cay.) ae 


a 


Q 
ct 


for the time (16 October 1973) under study. The net lateral 
flux-divergence term ee has not been assessed here, but pre- 
sumably its mean effect would show up as cold advection 
across the equator, and would result in a reduced value of 
S,. 

In Fig. 13, the dashed portion of the (E+H, ) and oe, oo) 
curves are displayed to show the effects of varying the K* 
value in the calculation of latent heat transport (Eq. 5-1), 
as discussed in Section V. In the dashed portion. the K* 
Twiewore 2 or lO- sem Sees’ was tested for latitudes 0°, 5°, 
Memo wavith aAstransitional value of 1.5 x 10° cm* sec”? 
used at 15°N. At all points on the solid (E+H,, ) line K* was 
1 x 10° cm* sec7™!. This smaller K* value seemed to provide 
more realistic values of (E+H,) in mid-latitudes but to some- 
what unrealistically small values in low latitudes. 

A possible explanation for the small (E+H,) values at 
low latitudes is the difficulty in correctly determining the 
temperature and humidity gradients between the surface and 


Uiewtop,of the constant flux layer in the data-sparse 
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equatorial areas. In these tropical latitudes some time- 
and space-smoothing of the few available radiosondes could 
have been introduced into the original FNWC map analyses, 
which would tend to smooth out vertical gradients as against 
those to be expected with observed radiosondes on 16 October 
1973. Hence it was considered justifiable to revert to 
Marbpaleto (1074) sctieeested K* = 2 x 10° cm* sec™’ in low 
latitudes, 0°-10°N, where the resultant dashed (E+H,, ) curve 
of Fig. 13 qualitatively appears more reasonable. For ex- 
ample, the peak in (E+H,, ) at 15°N appears to be due to rela- 
tively strong evaporative and sensible heat transfers 
associated with the Trade Wind regime, while the higher 
latitude peak (35°-50°N) in Fig. 13 is due to the trend to- 
wards winter maximum turbulent-heat transfers associated 
with cool continental air flowing over warm ocean currents. 
No attempt to select a final K*¥ value is suggested from 
these results. If the same turbulent-transfer model can be 
tested over a large enough set of seasonal samples, it should 
be possible to assess the appropriate K*-value which yields 


Pie DeSteabproximacilon to the atmospheric heat balance. 


D. SURFACE HEATING RATES 

ae 14 depicts the zonally-averaged surface heating bud- 
get as a function of latitude resulting from this study. 
The parameters include the net surface radiative warming (R), 
the surface heat loss due to turbulent transfer, -(EtH,), 
and the net warming or cooling at ae Surrace . Ce?) 


These parameters are related by the following expression 
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Figure 14. Surface heat budget disposition. The net warming 


(cooling) rate at the earth's surface no is given by the 
sum of the net radiative surface warming rate R and the loss 


rate due to turbulent heat transfér, -(E+H). 
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Table V. Comparisons of net radiative warming rate at the 
Surface (R), the heat loss rate at the surface due to turbu- 
lent transport -~CETH, ), and the net warming or cooling at the 
earth's surface (Qt) for 16 October 1973 (denoted by sub- 


script "1"). For comparison, annual mean values reported by 
Malkus (1962) are listed at selected latitudes (denoted by 
Sipseuupue 2.) l values are in ly min~?. 
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Qo * S = R - (EH) (7-4) 


where S is the oceanic heat-storage term. Again it was not 
possible here to separate the Pipa ions. oa and S. The 
Square points and dashed-in portion of the -(E+H,) and Q.. 
eUmves depict the Cftect Of varying the K* value in Eq. (5-1) 
momCisecussedeimeoces VileC, fhe overall distribution in the 
net oceanic heating rate cr) ) is the same for both choices 
of K*. However this rate is somewhat reduced for the case 

of the larger K* values suggested as applicable for the 
tropics. 

Table V shows quantitatively the difference in the param- 
eters depicted in Fig. 14 as found by this study for 16 
October i973 and the climatological annual values reported 
by Malkus (1962). Note the smaller mean seasonal radiation- 
al warming rate (R) and the effect of the smaller mean heat 
loss by turbulent transfer in giving a slightly larger mean 


warming of the surface. 
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VIII. CONCLUSIONS 


A heat budget model for inclusion in the FNWC numerical 
weather prediction scheme was evaluated uSing data along 
four meridional lines over the Pacific and Atlantic Oceans. 
These data soundings were used to compute a Space-time aver- 
age of the radiational and turbulent heating parameters for 
one day over an ocean environment. This space-time average 
agreed reasonably well with values derived from other stud- 
ies, including satellite measurements. There were some ob- 
vious differences which could be ascribed to the seasonality 
of the data used in this study comprising only a one day 
Space-time average while the comparison data was averaged 
over years to arrive at climatological results. 

Other differences could be ascribed to the two-cloud 
layer parameterization for cloud amounts, the treatment of 
cloud reflectivity and absorptivity for solar wave lengths, 
and the parameterization of sensible and latent heat trans- 
fers. At this stage in the model testing, it would seem 
that the radiational estimates are still somewhat question- 
able due to the adoption of the Smagorinsky-type cloud 
amounts. However the regular variation of the radiative 
terms from gridpoint to gridpoint tends to support the ap- 
parent realism of the radiative model. The turbulent trans- 
fer heating rates of the model seem to be subject to a 
greater degree of eee hae random error due to uncertainty 


of the vertical gradients involved near the surface. 
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To complete the comparisons of the radiational aspects 
of this model with cited climatological data (see Sec. VII) 
and with satellite data, further testing should Hemcdone at 
different seasons of the year and possibly over other oceanic 
areas. Thus with a full year of representative data, refine- 
ments of the cloud-cover parameterization and in the cloud 
absorption and reflection values used in the model may be 
assessed. 

Miialivemmilicsecbeatest pridpoint to gridpoint variability 
SecurbeasiiecOonnececion With the computation of evaporative 
and sensible heat parameters. It is perhaps here that the 
moisture and heating source terms used in the FNWC predic- 


tion model might best be improved. 
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